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METHOD USED FOR THE DETERMINATION OF THE 


EXPLORER XXVI SPIN AXIS POSITION 

P. A. Bracken 
L. R. Davis 
R. W. Janetzke 


ABSTRACT 

Two series of least squares fits were performed to determine the Explorer 
XXVI spin axis position in celestial coordinates throughout orbits 1 to 533, The 
data input to the calculations consisted of the sun’s position, Jenseiy-Cain 
magnetic field positions, the sun- spin axis angle, values of the magnetic field- 
spin axis angle and values of the angle between the spin-sun and the spin-field 
planes. The latter three angles were derived from the on-board Solar Aspect 
detector and from the roll modulated output curve of the Ion and Electron par- 
ticle experiment. 

The results of the fits are presented as a table listing a spin axis right 
ascension-declination pair along with the uncertainties in these angles for each 
orbit through orbit 533. 
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A METHOD USED FOR THE DETERMINATION OF THE 
EXPLORER XXVI SPIN AXIS POSITION 


INTRODUCTION 

This report describes a method used to find the Explorer XXVI spin axis 
position in celestial coordinates as a function of time through orbit 533, 

The determination was based on knowing the sun’s position, the magnetic field 
position (Jensen- Cain field, Reference [4]) and values of the angles: 

1. A. -the angle between the spin axis and sun positions, 

2. u.-the angle between the spin axis and magnetic field positions, and 

3. 0- the angle between the spin axis-sun and the spin axis-field planes 
measured in the direction of spacecraft rotation. 

Data values for the angle k were gotten directly from the Explorer XXVI 
Digital Solar Aspect Sensor (as described in Reference [1] ). Data values of the 
angles <x and 0 were determined from the roll modulated output curve of the 
Explorer XXVI Ion- Electron particle experiment (as described in Reference [8]). 

A non-linear least squares fitting program GLSWS was used in two separate 
series of fits to obtain the final spin axis coordinates. During the first series 
of fits, the sun and magnetic field positions along with values of k, a and 0 from 
one orbit were used to obtain a determination of the spin axis position. Because 
of the data selection criteria used in the first series of fits, from four to twenty 
separate determinations were possible for each orbit through orbit 540, 

A second series of GLSWS fits were performed fitting parabolas to the re- 
sults of the first series of GLSWS determinations versus orbit number. The 
parabola fits were performed on fifteen orbit batches. The final spin axis co- 
ordinates of the mid orbit of a fifteen orbit batch were taken as the value of the 
fitted parabola at that mid orbit. The fifteen orbit batches were slid up one 
orbit at a time allowing final spin axis determinations for orbits 7 through 533. 


1 
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1} CHAPTER I 

I PRELIMINARY CONCEPTS 


A. THE BASIC PROBLEM 


The objective of the study described in this report was to determine the 
Explorer XXVI (EPE-D) spin axis orientation as a function of time in inertial 
coordinates. The determination was performed on orbits 1 through 533, a period 

during which Explorer XXVI was spin stabilized. * 

i 

ij The data used in the determination consisted of: 

- 

'f. 

i] X. 

U 
.?! 

t .. 2 , 

) '■ 

j( 

3 " 3. 

I 


The angles k, a and 0 are shown in Figure 1 for one possible sun, magnetic 
field, spin axis orientation. 


measurements of the spin axis- sun angle, k t provided by the on-board 
solar aspect sensor, 

values of the angle (0) from the sun-spin axis plane to the magnetic field- 
spin axis plane measured in the direction of spacecraft rotation, 

values of the angle (a) between the magnetic field and spin axis positions, 
and 

the statistical uncertainties cr K , cr, and cr a in k , i p and a. 


U 



The magnetic field position in celestial coordinates was that given by.the 
Jensen and Cain Magnetic Field Coefficients for 1960 (see Reference [4]). 

The sun's position in celestial coordinates was derived by using the com- 
puter program HE LOS as described in Reference [6] . 

Values of the a and </> angles were derived from solar aspect data and from 
the Ion and Electron detector's roll modulated output curve. For a description 
of the I&E detector see Reference [8], 



As is illustrated in Figure 2, k and a fix the spin axis position to the inter- 
sections of circles of constant a and of constant k on the celestial sphere. If in 


♦Explorer XXVI was launched on December 21, 1964. Orbits 1-533 cover the approximate time 
period from day 356.3927 (December 21, 1964) to day 524.9595 (June 7, 1965). 
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addition, tho angle # is known then ambiguity Is eliminated and the true spin 
axis pv f'ition (l, e. either SA t or SA 2 in Figure 2) can be determined, (The di- 
rection of rotation of the spacecraft was known.) 


In general, knowledge of any two of the angles A, a and </> will yield two pos- 
sible spin axis positions and knowledge of all three of the angles will resolve 
the ambiguity. In practice, there are eases when knowledge of only two of A, a 
or it is sufficient for a spin axis determination because } 

1. the spin axis location moves slowly with time and its position is known 
well enough to resolve ambiguity, and 

2. a number of a and 4* values can be used over periods of up to one hour 
during which time the field direction will change significantly thus help- 
ing to resolve the spin axis ambiguity, 


B. RELATIONSHIP OF 0 AND a TO THE I&E DETECTOR 
OUTPUT CURVE 


Charged particles having small pitch angles* are lost into the atmosphere. 
For this reason, the I&E detector will measure a minimum in particle intensity 
when looking along a magnetic field line. Similarly, a maximum in particle in- 
tensity occurs when the I&E detector measures particles having 90° pitch angles- 
that is, when the I&E detector is pointing in a direction 90° away from the field 
line. This response was found to be empirically true for the I&E detector at 
selected times and was assumed to be true for all times. 

On Explorer XXVI, the I&E detector was mounted at a 45° angle to the 
spacecraft spin axis so that a cone of half angle 45° was swept out about the spin 
axis by the I&E detector during each roll. By making use of the solar aspect 
data provided on Explorer XXVI,. the output of the I&E detector can be plotted 
against the fraction of a roll after see sun at which each I&E measurement was 
taken. f 

A direct reading from such a roll modulated output curve of the abscissa 
corresponding to a minimum in particle intensity will yield either i> or i p & 180°, 
Figure 3 illustrates the point for one possible roll modulation pattern. In Fig- 
ure 3, the abscissa % x represents the fraction of a roll after see sun at which 
the I&E detector was looking in either the plus field or the minus field directions 


•The pitch angle is the angle between the particle velocity vector and the magnetic field vector 
at a point in space. 
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(Similarly lor the abscissa x 2 ), Together and x 2 correspond to one pair ol 
values from the soft#, $ + 180% # - 180°}* 

The relationship of a to the I&E roll modulation curve can be seen by ex- 
amining Figure 4. The Taw of Cosines applied to the spherical triangle in Fig- 
ure 4 with arcs a, 90° and 45° yields? 


cos 90® s cos 45° * cos a * sin 45° * sin a. • cos (180° ~ f/2) (1.1) 


clLuj 

mm 


where f is the fraction of a roll between maximum particle intensity readings by 
the I&E detector during one roll. 


Solving Equation 1. 1 for a gives 


a 


“ 90° - a ret an (cos (f/2)) 


( 1 . 2 ) 


Thus it is seen that a is a function of the fraction of a roll between maximum 
particle Intensities as measured by the I&E detector. 


As with 'A# the determination of a from the I&E output curve can yield either 
a or 180° - a. Thus the a, \p data used in the spin axis determinations are am- 
biguous in that it is not known whether they are measured with respect to the 
positive or negative field direction at any given time. 


Examples of possible I&E detector output curves are given in Figures 5 
and 6. 


Further information on the actual method used to calculate a and $ from the 
I&E detector data is given in Reference [3] . Suffice it to say that in practice 
V> can always be found from the I&E data but that a cannot be calculated outside 
the range 60° < a < 120°. When a is less than 60° or greater than 120% maximum 
intensity occurrences are either too flat or otherwise too ill defined on the I&E 
roll modulation curve to be used in a determination of a. 


C. DETERMINATION OF THE SPIN AXIS-SUN ANGLE 


The solar aspect sensor on board Explorer XXVI provided digitized meas- 
urements of the sun-spin axis angle, k . The solar aspect detector, with a fan 
field of vision 180 degrees in length, was mounted so that the fan swept out the 
entire sky as the satellite rotated. 















r 

? i 

i 

5 The nominal uncertainty in the sun-spin angle measurements was £L. ii * 

since the solar aspect detector sampled the 180° fan in 63 approximately equal 
intervals. However when the X -reading changes from one of the 63 sectors to 
another, the angle is known more accurately— to approximately * 0 . 3 degrees. 
Also when the motion of the spin axis is regular with respect to the sun, the 
angle X may be interpolated to this higher accuracy. 

For further information on the solar aspect detector see References Cl] 
and [ 2 ], 



CHAPTER ll 


METHOD FOR SPIN AXIS DETERMINATION 


A. THE MATHEMATICAL PROCEDURE FOLLOWED 

. Denote the right ascension and declination of the sun, magnetic field and 

spin axis by SRA, SDC, FRA, FDC, ARA and ADC respectively. Using this notation, 
the angles X. , a and 0 can be written as 

f X s f(ARA, ADC, SRA, SDC) 

• y •' 

y - a = g(ARA, ADC, SRA, SDC, FRA, FDC) 

0 “ h(ARA, ADC, SRA, SDC, FRA, FDC) 

I (2.1) 

, f _ V - 

/•••' i 

V 5 where f, g and h are unknown functions. 

0; For a given time period the following data are known: 

:r -a 

__ a. the sun- spin axis angle A 

B\ 

X \ b. the celestial coordinates SRA, SDC of the sun 

ll c. n values of 0 as P and the associated field coordinates FRA it FDC. . 

M 

%i d. m values of a as A. and the associated field coordinates values FRAj , 

1! FDC., and 

y:.h J 

1 ‘ ' ' " . . 

W‘ e. the statistical uncertainties cy> .}" =1 and { cr Aj} m -i * or die angles 

yj in (a), (c) and (d) respectively. 11 3 3 >■> 

The problem was to find those values of the spin axis coordinates ARA and 
ADC which when used in Equations 2. 1 along with the proper field and sun co- 
ordinates would yield a minimum value for the expression: 


-■ • !■*. 

I I 



where X, \jt i and a. are the calculated values yielded from Equations 2.1 and 
A, 'l' . , Aj are the measured data values. 

Given the data as listed in (a) -(e) above, a modified version of W. E. Daniels 
General Least Squares (GLSWS) program was used to solve for the spin axis co- 
ordinates during given time periods. A description of GLSWS is given in Ref- 
erence [7] . 

For ease in using GLSWS, Equations 2. 1 were rewritten in one equation as 


y = F(ARA, ADC, SRA, SDC, FRA, FDC, I) (2.3) 


where the independent variable I is 1, 2 or 3 for the function F to yield values 
ofX, a or 0 respectively. 

In conjunction with the notation used in Reference [7] , Equation 2. 3 was 
written as 


y “ f(p1, P2, Xj, X 2 , X 3 , X 4 , x s ) 


(2.4) 


where Pi andP2 are the spin axis right ascension and declination and x 1 through 
x s are SRA, SDC, FRA, FDC and I respectively as in Equation 2. 3. 


The program GLSWS uses an iterative algorithm to accomplish the required 
fit. To initialize the algorithm, it was required that the following values be fed 
to GLSWS: 


a. guess values of ARA and ADC as Pl 0 and P2 0 

b. k data values Y. and (x t . , • • * , X s . ) 


c. k calculated values y. (derived from Equation 2.4) using the values 
(x u , ••• , x 51 ) andPl 0 ,P2 0 


d. the partial derivatives 
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(JJL) 

\dP2j 


Pl 0 ,P2 0 ,X n ,-'V,X 5: 


(2.5) 


for i « 1, • • • , k. 

In return, GLSWS yielded new values of Pi and P2 which minimized the 
expression 




k ~ 2 


( 2 . 6 ) 


Since the form of the function F in Equation 2.4 was quite long and compli- 
cated, a practical method of calculating the values in (c) and (e) above was 
evolved. This method consisted in utilizing the known celestial coordinates of 
the sun and field positions together with the guess coordinates of the spin axis 
position to obtain spherical triangle solutions yielding calculated values y . for 
X, a and For example, in Figure 7 the codeclination of the sun and spin axis 
along with the difference in their right ascensions provides enough information 
to solve the given spherical triangle for k. Similarly, spherical triangle solu- 
tions may be obtained for a and \Jj. In practice, the program CAAAB as described 
in Reference [5] was used to solve the required spherical triangles. Using the 
above procedure, spherical triangle solutions were found for the k required 
calculated values: 


F(Pl 0 ,P2 0 ,X li , X s .) . 


(2.7) 


In addition, the guess values Pl 0 andP2 0 were incremented by an amount A* and 
these incremented values were used in spherical triangle solutions to obtain 
the calculated values: 


y' = F(P1 0 +A, P2 0 ,X U , ‘ * t X Si ) 


F(F1 0 ,P2 0 +A, X n , •••,X S1 _) 


( 2 . 8 ) 


*In our calculations A was taken to be one degree. 
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Finally the required partial derivatives were calculated from: 


f dF_ \ 
\dPl) 


Vi "Vi 


Pl 0 ’ P2 0 ' X li' — • x s i 


A 


(JL) 

\dP2/ 


y/"Vi 


P»0' P2 0 ' X 1 i * * * * ' x s i 


A 


( 2 . 9 ) 


At the beginning of the n + 1 st iteration, the values y i and the required 
partial derivatives were calculated from Equations 2.7, 2. 8 and 2. 9 using the 
values PI „ andP2„ of the spin axis coordinates returned by GLSWS at the end of 
the n iteration. This procedure was allowed to continue until successive 
iterates were close enough (difference <10“ 2 degrees) or until the number of 
iterations exceeded a specified value. 


B. THE METHOD AND INPUT DATA USED 


In practice, the mathematical procedure and the da;a input to GLSWS were 
modified slightly to force GLSWS to perform a weighted least squares fit utilizing 
the known uncertainties cr a and cr 


To accomplish the weighted fit, the basic equation (Equation 2. 7) was utilized 
to obtain a new set of calculated values as: 


y i - Y i 


z . 

i 


cr. 


Y , 


(2.10) 


where Y . , cr . are the measured data values and their uncertainties and where 

V * i 

y i is given by Equation 2.7 for i - 1, . . . , k. 


Redefining the basic expression as in Equation 2. 10 caused GLSWS to per- 
form the weighted fit while minimizing the expression: 


RMS 


/EE 

~ y k -2 


-0.0) 

2 
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I { 
1 4 


Y- \- 

r { 
? : 

5 • v 


{ i 


i •, i 

\ 4 


l ' V 


1 1 
! v 


or 


» 


RMS 




k-2 


( 2 . 11 ) 


Also the partial derivatives supplied to GLSWS became: 


and 






Pt n ,P2 n ,X u ,--,X sl 


( 2 . 12 ) 


(2. 13) 


where i « 1, » . . , k. 

hi summary, to initiate a fit GLSWS was supplied with: 

t 

a. k values of Y . , o* y and (x u , * * * , X 5 { ) 

b. guess values Pl 0 andP2 0 for PI andP2 

c. k values of z. calculated from Equation 2.10 and 

d. k partial derivatives dz ./&P1, dz L /dP2 as calculated from Equations 2. 12 
and 2.13 using (Pl 0 , P2 0 ) for (pi n , P2 n ). 

The data supplied to GLSWS for a spin axis determination were selected 
according to the following criteria: 

1. the data for each fit had to be from the same orbit 

2. the data for each fit had to extend over a time period of no longer than 
one hour 

. ‘11 > f 'd; : : " 
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3, the data used for a fit had to extend over an interval of no larger than 
. 5 earth radii in L value (ex. data from L «* 3. 0 to L ~ 3. 5) and 

4. the total number of points used for a fit had to be at least four but less 
than ninety-seven. 

For any given fit during which k data values were supplied to GLSWS, k was 
comprised of one value of k, n values of a and k - n “1 values of 0 where n 
varied from 0 to (k - l)/2. 

The uncertainties and cr, in a and 0 were fixed at three degrees and two 
degrees respectively. The data values for a and 0 were read in from magnetic 
tapes whose format is described in Reference [3] . 

The data values used for k and cr k are tabulated in Table X. The gaps in 
this table were filled in by linear interpolation. 
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CHAPTER HI 


CONCLUSION 
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A. PRESENTATION AND DISCUSSION OP THE RESULTS DERIVED 

As a result of the criteria applied during input data selection, there were 
as many as twenty separate spin axis determinations made by GLSWS for each 
orbit in the range from orbit 1 to orbit 540. 

A plot of the twenty or so spin axis determinations for each orbit from the 
first round of GLSWS fittings shewed considerable scatter in these results. To 
eliminate this scatter and to achieve the most reliable final result, an additional 
fit was performed on the spin axis coordinates determined from the first round 
of GLSWS fits. In this second GLSWS fit, parabolas were fitted to the spin axis 
coordinates from the first GLSWS determinations on fifteen orbit batches. The 
coordinates yielded from the fitted parabolas for the mid orbit of each fifteen 
orbit batch were taken as the final spin axis position for that orbit. Also the 
fifteen orbit batches were slid up one orbit at a time so that the coordinates for 
orbit n were derived from a fit to the GLSWS determinations for orbits n - 7 
to n + 7. 

In addition, only those original determinations (of the twenty or so possibil- 
ities) which resulted from fits having at least ninety percent as many a’s as 0's 
and which arose from data having L values between 1. 5 and 4. 5 were used in the 
final parabola fits. Consequently, only from three to ten of the twenty possible 
original determinations per orbit were used as input data for the fifteen orbit 
parabola fits. 

The results of the final parabola fits are listed in Table 2, Columns three 
through eight of Table 2 are expressed in degrees* Column two lists orbit start 
times in day plus fraction of the day where the following conventions hold: 

X. launch occurred on December 21, 1964 (day 356* of 1964) 

2. day 367. corresponds to January 1, 1965 (since 1964 was a leap year) 


3. day 524. corresponds to June 7, 1965. 


Columns three and four of Table 2 list the final values of the spin axis right 
ascension and declination respectively. Columns seven and eight are the final 
values of the parabola fit RMS yielded by GLSWS (as expressed in Equation 2. 6). 
Columns five and six are the uncertainties in columns three and four respectively 






and were determined from: 


^ra s ^ra ' r D n 

^DEC “ ®M®DEC ‘ M (3.1) 


Djj and D A j are the first elements of the inverse coefficient matrices from the 
final parabola fits in GLSWS. 

The coordinates listed in Table 2 for orbits one through seven are '’eyeball’' 
extrapolations obtained from a plot of the remainder of the table from orbit eight 
on. Consequently, the uncertainty and RMS values in columns five through eight 
for the first seven orbits are listed as "999, ", 

A plot of the results listed in Table 2 appears in Figure 8. 


B. LIMITATIONS OF THE METHOD USED 

In some orientations of the spin axis, sun, field positions one of the spin 
axis coordinates can be heavily dependent on one of X., a or 0. Thus for a good 
spin axis determination in such cases, values of that angle must be present in 
the input data fed to GLSWS during the first round of fits. For example, from 
orbit 400 on the Explorer XXVI spin axis declination was ill determined when 
only values of \ and 0 were used for a determination. For this reason, only 
those spin axis determinations calculated for time periods with X, a and 0 values 
present in the input data were used to obtain the final results, For the calcula- 
tions from orbit 400 on, best determinations were made when there were almost 
as many a's as 0's present in a block of data fed to GLSWS during the first round 
fits. 

It was also found that data associated with L values greater than 4. 5 some- 
times yielded erroneous spin axis determinations. For this reason the input 
data used were restricted to L shells in the range 1. 5 < L < 4.5. 

The erroneous spin axis determinations probably resulted for one or more 
of the following reasons: * . 

1. at certain times the effects of magnetic disturbances were still signifi- 
cant down to L shells of L - 4. 5, or 
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2, flat pitch angle distributions were encountered yielding inaccurate de- 
terminations of a and 4> from the I&E detector output, or 

3. count rates went down while statistical fluctuations became 
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APPENDIX A 
NOTATION 


A, A ~ the spin axis-sun angle, 

a, a - the spin axis-magnetic field angle, 

0, f - the angle between the spin- sun and the spin-field planes measured in 
the direction of spacecraft rotation, 

tr k , - the statistical uncertainty in A, A, 

er a , - the statistical uncertainty in a, A. 

0 \j f - the statistical uncertainty in <//, f , 

ARA - spin axis right ascension, 

ADC - spin axis declination. 

SRA - sun right ascension, 

SBC - sun declination. 

FRA - magnetic field right ascension. 

FDC - magnetic field declination, 

y . - a calculated value of A, a or 1 /'. 

Y L - a measured value of A, a or 0, 
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APPENDIX B 


TABLES 


Table 1 


Values of the Spin Axis-Sun Angle 
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Figure 5. I&E Roll Modulation Curves for Values of a up to 90° 
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Figure 7. Calculated Values for X. 
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Figure 8b. Explorer XXVI Spin Axis Coordinates 
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